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Summary 

The Five-Axis, Three-Magnetic-Bearing Dynamic Spin 
Rig, a significant advancement in the Dynamic Spin Rig 
(DSR), is used to perform vibration tests of turbomachinery 
blades and components under rotating and nonrotating 
conditions in a vacuum. The rig has as its critical compo- 
nents three magnetic bearings: two heteropolar radial ac- 
tive magnetic bearings and a magnetic thrust bearing. The 
bearing configuration allows full vertical rotor magnetic 
suspension along with a feed-forward control feature, 
which will enable the excitation of various natural blade 
modes in bladed disk test articles. The theoretical, me- 
chanical, electrical, and electronic aspects of the rig are 
discussed. Also presented are the forced-excitation results 
of a fully levitated, rotating and nonrotating, unbladed rotor 
and a fully levitated, rotating and nonrotating, bladed rotor 
in which a pair of blades was arranged 180° apart from each 
other. These tests include the “bounce” mode excitation of 
the rotor in which the rotor was excited at the blade natural 


frequency of 144 Hz. The rotor natural mode frequency of 
355 Hz was discerned from the plot of acceleration versus 
frequency. For nonrotating blades, a blade-tip excitation 
amplitude of approximately 100 g/A was achieved at the 
first-bending critical (-144 Hz) and at the first-torsional 
and second-bending blade modes. A blade-tip displace- 
ment of 70 mils was achieved at the first-bending critical by 
exciting the blades at a forced-excitation phase angle of 90° 
relative to the vertical plane containing the blades while 
simultaneously rotating the shaft at 3000 rpm. 


Symbols 

A cross section 
c constant defined in eq. (A9) 
B magnetic flux density 
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F net controlling force on rotor 
/ force on rotor due to one pole 
g acceleration due to gravity (9.8 m/s 2 ) 
H magnetic field 

i current in coil 

K control gain 

X stiffness 

k an rigid rotor equivalent spring constant 

lj e mean length of magnetic displacement 

m pn rigid rotor equivalent lumped mass 

n number of turns in coil 

5 air gap spacing 

Sfc linear velocity of shaft 

Sfa linear acceleration of shaft 

V volume 

W magnetic energy 

<E> magnetic flux 

(I permeability 

q 0 permeability of free space 

CO excitation frequency 

Subscripts: 

A a through air gap cross section 
A fe through iron cross section 
a air gap between stator and rotor 


b bias 

c control 

d derivative 

de deviation from nominal gap distance between stator 
and rotor 

eq equivalent 

ex excitation 

fe iron 

i current 

/ linearized 

p proportional 

r relative 

s position 

0 nominal 

Introduction 

Vibration is a bane in turbomachinery and is responsible 
for more than 40 percent of the blade failures in gas turbines 
(refs. 1 to 3). Blade vibration in jet engines is responsible 
for noise and blade fatigue. Fatigue poses a safety risk; 
therefore, it is important to devise mechanisms and tech- 
niques for studying vibration and its effect on 
turbomachinery components. The NASA Glenn Research 
Center has been working on structural vibration and con- 
trols in turbomachinery and has built the Dynamic Spin Rig 
(DSR, ref. 4) to assist in this research (fig. 1). 

The DSR employed two voice-coil-type linear electro- 
magnetic shakers to apply oscillatory axial forces to the 
rotor and test blades through a thrust bearing. A subsequent 
upgrade of the DSR (ref. 5) incorporated a single heteropo- 
lar radial active magnetic bearing, which allows magnetic 
excitation and two-axis radial suspension of the rotor. The 
upgraded DSR (fig. 2) supported a vertical rotor by 
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Figure 1 . — Dynamic Spin Rig. 


employing a conventional radial ball bearing at the upper 
end and a radial magnetic bearing at the lower end. 

The rig shown in figure 2 was used to perform vibration 
tests of turbomachinery blades and components under 
spinning conditions in a vacuum. However, the ball bearing 
imposed limitations on the maximum rotational speed of 
the rotor (18 000 rpm) and on the test duration at lower 
speeds due to frictional heating. These operational limita- 
tions adversely affected the nature and type of research that 
could be conducted on test articles. Indeed, it was the desire 
to surpass the rotational speed and operational duration 
ceilings and to have the enhanced excitation capabilities 
that gave rise (ref. 5) to the late 1990’s notion of total 
magnetic suspension for the DSR. 

The Five-Axis, Three-Magnetic-Bearing Dynamic Spin 
Rig is the progeny of the total rotor suspension goal and 
represents a significant advancement beyond the capabili- 
ties of the previous rigs. This advanced rig (figs. 3 and 4) 
is discussed in detail in this report and consists of a 
magnetic thrust bearing, a magnetic radial upper bearing, 
and a magnetic radial lower bearing. Both radial bearings 
are heteropolar. Simultaneously energizing the bearings 
under closed-loop control (refs. 6 and 7) effects full levita- 
tion of the rotor, making it possible to achieve higher 
rotational speeds (in excess of 18 000 rpm). The upper 



Figure 2. — Upgraded Dynamic Spin Rig showing conven- 
tional upper bearing and lower radial magnetic bearing. 


rotational speeds are now limited mainly by containment 
issues and stresses in the test articles. The higher rotational 
speed is a critical element in high-cycle-fatigue research 
projects pertaining to the damping and mistuning of bladed 
disks. This rig design should also allow higher excitation 
amplitudes (by virtue of full rotor suspension, which permits 
larger rotor translation and tilt displacements) than were 
achievable with the previous rig configuration. 

This report presents the basic theory and design parameters 
of the new rig and also discusses the excitation data gleaned 
from tests involving (a) an unbladed, nonrotating and rotating 
(10 000 rpm) rotor and (b) a bladed rotor that was excited in 
the nonrotating manner and then subsequently spun at 
3000 rpm. The next section discusses the requisite com- 
puter hardware specification for successful rig operation. 
The reader is directed to appendixes A and B for additional 
information related to the theoretical aspect of the rig. 
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Figure 3. — Engineering schematic of five-axis ensemble. (All dimensions are in inches.) 


Computer Hardware 

Experience has shown that robust control of the five- 
axis-rig rotor at all operating speeds requires a control code 
loop time at or below 100 |is. To achieve this level of 
performance, a 400-MHz Pentium II processor was used to 
effect a code loop time of 50 |is. There are 8 input board 
channels (5 proximeter sensors, one 24-per-revolution, one 
1-per-revolution, 1 external signal import) and 11 output 
board channels. Of these, four are used for the lower 
bearing, four for the upper bearing, and two for the thrust 
bearing. The remaining channel is used for outputting the 
optional excitation signal (generated within the code) to an 
oscilloscope. The bearings and bearing components are 
discussed in the next section. Reference 6 provides addi- 
tional hardware information and a detailed explanation of 
the 1-per-revolution signal logic. 


Bearing Components and Rig Support 
Structure 

The five-axis rig contains three magnetic bearings: a 
thrust bearing, a radial upper bearing, and a radial lower 
bearing. The thrust bearing has 2 thrust coils, each of which 
has 225 turns of 19-gauge copper wire. Each radial stator 
has 4 coils, each having 156 turns of 2- by 24-gauge copper 
wire. Figure 5 shows the upper and lower radial stators; 
figure 6, the axial view of a radial magnetic-bearing stator 
lamination showing four three-sectional pole elements; 
figure 7, the actual rotor hanging from the rig support 
structure; figure 8, the top portion of the rotor where the 
thrust bearing is affixed. Figure 9 illustrates the thrust plate 
and the thrust coils, and figure 10 depicts the fully as- 
sembled rig supported by the test cell service tripod. In 
high-speed experiments, the rig is operated in the vacuum 
chamber shown in figure 1 1 . 
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Figure 4. — Five-axis rig in its current configuration. 



Figure 7. — Rotor hanging from rig support structure. 
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Figure 8. — Top view of rotor without thrust bearing 
assembly. 



Figure 9. — Thrust plate and thrust coils. 


It should be noted that the multisectional design of each 
of the four poles in the radial stator (fig. 6) is dispensable. 
These radial laminations were adapted from a motor project 
as a cost-cutting measure. The unusual bore geometry in 
the thrust plate (known as a three-point polygon, fig. 9) 
obviates the need for a shrink fit onto the rotor, thus 
minimizing radial and tangential stresses between the plate 
and rotor. The top portion of the rotor, having a similar 
geometry but a slightly smaller dimension, fits snugly into 
the bore. 

A commercially available power amplifier (coupled 
with the aforementioned computer controller) is used to 
energize and control each coil in the bearings discussed 
above. Each coil is driven by a pulse-width-modulated 
(PWM) amplifier that can output up to 170 V at 15 A 
continuously and 30 A for short durations. The amplifiers 
and the tuning procedure are discussed next. 



Figure 10. — Fully assembled rig mounted on test cell 
service tripod. 



Figure 1 1 . — Vacuum chamber. 
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Power Amplifier Tuning 


10.00 


For stability and maximum bandwidth, each amplifier- 
coil pair has to be matched. Prior to this task, tests were 
performed to detect wire-to-wire or coil-to-stator shorts. 
These tests include resistance and inductance measure- 
ments, continuity, and megger wire insulation tests. A low- 
amplitude dc signal was applied to the actuator coils, and 
the flux polarity was determined by using a thin Hall-effect 
probe inserted in the air gap between the center of the pole 
and the rotor. Knowledge of the flux direction is necessary 
for the apposite design of the control code voltage output. 
These and some tuning tests were conducted while main- 
taining a uniform air gap (using the mechanical clamps of 
the rig) between the stator poles and the rotor. 

NASA Glenn researchers modified the PWM amplifiers 
by replacing a fixed resistor in the R-C error amplifier 
compensation circuit with a 0- to 2-MQ variable resistor. A 
swept-sine signal was inputted to the amplifiers and the 
ratio of the current sensor output to the swept-sine voltage 
input was monitored in the frequency and phase domain. 
The 2-MQ potentiometer was adjusted to optimize ampli- 
fier bandwidth while maintaining stability (see figs . 1 2 and 
1 3 for the resulting transfer function plots of the coils in the 
bearings). It should be noted that the two coils in the thrust 
bearing are nearly identical and the eight coils in the radial 
bearings (four in each bearing) are nearly identical; hence, 
only two typical transfer function curves are presented: one 
for the thrust bearing and one for the radial bearings. 

Figure 12(a) shows that the magnitude of the power 
amplifier gain (current out divided by voltage in) for the 
thrust coil is approximately flat for input excitation signals 
ranging up to about 2000 Hz. The accompanying phase- 
shift curve (fig. 12(b)) indicates that at 1000 Hz, the output 
current is phase shifted -45° relative to the input signal. For 
the radial bearings (fig. 13), the gain is fiat to about 1000 Hz 
and the output current is phase shifted -90° relative to the 
input signal. 

Between each amplifier and its associated bearing coil is 
a commercially available current sensor (fig. 1 4). Using the 
current sensor output and the four trim pots on the PWM 
amplifier, each amplifier was tuned to its associated coil. 
Direct-current offsets were eliminated and the amplifier dc 
gains were set to 1 A/V. A low-amplitude, 100-Hz sine 
wave was inputted to the PWM and the current sensor 
output monitored in the time domain to verify phase, 
amplitude gain, offset, and noise content. After fine adjust- 
ments were completed, other input sine wave signals at 
various frequencies and amplitudes were checked. 
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Figure 12. — Transfer function for thrust bearing, (a) Gain, 
(b) Phase. 


c 

0 

O) 

L_ * 
0 ! 
It ' 

Q. 

E 

< 


10.00 I— 


1.00 


0.10 



0.01 


180 

O) 

0 

■o 

90 

<D 

(A 

0 

-C 

CL 

0 

-90 


-180 


(a) 


10 


(b) 


10 


Source: +1 .50 -V bias 

0.25 -V swept sine 


100 1000 10 000 
Frequency, Hz 


100 1000 10 000 
Frequency, Hz 


Figure 13. — Transfer function for radial bearings, (a) Gain, 
(b) Phase. 


For the effective application of any excitation forces to 
the suspended rotor, the levitation of the rotor must be 
robust. It is therefore imperative that the stiffness-damping 
parameters be suitably adjusted. These parameters are 
discussed in the following section. 
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Power 

amplifier 

— Output 
filter 

— Current 
sensor 


Figure 14. — Bank of power amplifiers with their associ- 
ated output filters and current sensors. 


Stiffness-Damping-Parameter Stability 
Map 

Robust levitation of the rotor requires the appropriate 
setting of the digital (proportional-derivative) controller 
stiffness and damping values (ref. 6). For this report, 
however, the latest version of the code discussed in refer- 
ence 6 was employed to obtain the nonrotating stability 
limit curves for the thrust, upper, and lower bearings 
(figs. 1 5(a) to (c)). It is evident from these curves that for a 
given stiffness value, the corresponding damping values 
should reside between the extreme limits defined by the 
stiffness-damping curve envelope. For the current rig con- 
figuration, a stiffness gain of 1.20 digital units and a 
damping gain of 35.0 digital units were used for the thrust 
bearing. The upper bearing stiffness-damping gains were 
1.7 and 35.0 digital units and for the lower radial bearing, 
stiffness-damping gains of 1 .7 and 25.0 digital units. It was 
observed that these values were also appropriate for rota- 
tional speeds up to 10 000 rpm (the highest speed at which 
the rig was tested without blades). These digital units are 
specific to the controller code (ref. 6 and its latest version) 
and are the analogs of the stiffness and damping used in 
classical spring mass-damper dynamics. 



Stiffness, digital units 




Figure 15. — Static stability limit, (a) Thrust bearing, 
(b) Upper bearing, (c) Lower bearing. 


A number of novel features and refinements were imple- 
mented in this latest version of the controller code. These 
features include logic for (a) thrust bearing soft touchdown, 
(b) thrust plate demagnetization, (c) automatic angular 
sweep of the forced-excitation vector, and (d) adjustment 
of the tilt pivot point along the rotor. These features are 
discussed in detail in the following section. 
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New Software Features 

To avoid potential damage to the rotor and the axial 
backups, a soft delevitation scheme was incorporated in the 
controller. F igure 1 6 shows the controlled rampdown of the 
thrust plate position from the axial levitation point to the 
touchdown surface. This rampdown is achieved by com- 
manding a linearly decreasing levitation center by adjust- 
ing the offset term added to the force equation signal until 
touchdown occurs. Figure 1 7 shows the current level in the 
upper axial magnetic bearing during this controlled 
rampdown. Note the slight increase in the current value 
from 0.5 to 1.3 sec, which indicates the control algorithm 
is working to overcome the increased magnetic reluctance 
in the upper gap as it attempts to preserve levitation during 
the rampdown. 

Periodically, axial thrust disk remnant magnetization 
hinders axial relevitation using baseline controller gains. 
This situation usually requires the modification of the 
control bias in addition to the controller feedback stiffness 
to pull the thrust disk away from the top or bottom backup 
surface. To eliminate this problem, a new demagnetization 
scheme was implemented. In figure 17, the demagnetiza- 
tion current duration ranges from 1 .9 to 4.0 sec. At 1 .9 sec, 
the rotor is sitting just above the lower backup surface when 
the controller forces the current into an exponentially 
decaying cosine oscillation with a frequency of 5 Hz and no 
dc offset. This full-reversal-current signal, which decays in 
about 2 sec, significantly reduces the remnant field in the 
disk. In addition, the controller keeps track of the direction 
of the current in both axial stators during levitation. Upon 
delevitation, the controller reverses the direction of the 
current for the next operation. This logic helps counter any 
remaining remnant field not suppressed initially using the 
decay technique. 

Another novel feature implemented in the code is the 
user-defined and user-initiated automatic angular sweep of 
the forced-excitation signal used to excite the blades. The 
excitation angular sweep can range from 0° to 360° relative 
to the vertical plane containing the test blades as they rotate. 
However, an angular sweep from 0° to 1 80° is typical. The 
implementation of this option was motivated by the obser- 
vation of a diminution in blade response amplitude at 
speeds above 7500 rpm as measured by stationary fiber 
optical displacement transducers. For example, the excita- 
tion phase angle at maximum blade response had shifted 
from the desired 90° to a new value of 95 ° . The authors infer 
that this angular change is directly related to the phase shifts 
in the currents depicted in figures 12 and 13. 



Time, sec 


Figure 16. — Soft delevitation of Dynamic Spin Rig 
rotor. 



0 12 3 4 


Time, sec 

Figure 17. — Current in upper axial stator during 
delevitation and demagnetization. 

In an actual test (at rotor speeds at or above 7500 rpm), 
an excitation angular sweep is initiated over a set time 
interval (10 sec for example) and simultaneously the blade 
response (obtained with the aid of stationary fiber optical 
displacement transducers) is measured over the selected 
time interval. The angle at maximum resonance response is 
noted (at the desired rotor speed) and then is set for use in 
subsequent blade excitation and/or damping tests at that 
rotor speed. 

For flat-plate blades oriented in a vertical plane, the 
ability to excite torsional modes in the blades is a desirable 
feature for the rig. Torsional mode is best excited by rotor 
tilt motion. The tilt pivot point along the rotor may be 
adjusted upwards or downwards through the center of the 
hub or disk, allowing the operator to optimize the pivot 
location for maximum torsional effect. At the time of this 
writing, this option had not been exercised for the test 
blades. The excitation tests were conducted in the bounce 
mode only, the results of which are discussed extensively 
in the following section. The latest version of the code is 
available online at https://technology.grc.nasa.gov/software . 
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Demonstration of Excitation Capability 


As mentioned previously, the five-axis rig represents a 
significant advancement in the DSR because the rotor 
speed is not limited to 1 8 000 rpm and the rig operator now 
has the ability to excite (at various frequencies) rotor rigid 
body motions (bounce and tilt) for a more effective and 
thorough investigation of blade excitation and damping. 
The rotor excitation is achieved by applying an excitation 
signal concurrently to the upper and lower radial bearings. 
At each bearing, excitation signals are applied simulta- 
neously to the x- and y-axes. This simultaneous excitation 
produces a resultant force vector with a magnitude and an 
angular orientation. The direction and amplitude of this 
force vector can be fixed in a nonrotating laboratory frame 
of reference or made to rotate with the test article by 
enabling the synchronous-rotating-force feature in the con- 
trol code (ref. 6). A tiny mirror attached to the shaft reflects 
a beam of laser light once every rotation of the shaft. A 
sensor then converts the reflected light pulses to electrical 
pulses. The pulses are sent to an input channel on the A/D 
board and are used to trigger the synchronous-rotating- 
force logic block in the code, which then synchronizes the 
force vector with the rotating shaft. 

The excitation capability discussed above was demon- 
strated in the following set of experiments. In these tests, 
the rotor and a two-blade test article were instrumented 
with three Kistler accelerometers (100 mV/g) to character- 
ize the nonrotating excitation potential of the five-axis 
magnetic suspension system. Figure 1 8 shows two of these 
accelerometers: one mounted to the blade and one mounted 
on the blade-base support structure. The third accelerom- 
eter (not shown) was mounted perpendicular to the blade 
row at the bottom of the rotor. 

To characterize how much excitation capability exists, a 
determination of the amount of acceleration and displace- 
ment that can be provided at the base of the blade is 
appropriate. All data to follow were recorded using an HP 
3665A Spectrum Analyzer in the configuration shown in 
figure 1 8. Figure 1 9 presents curves of blade -base accelera- 
tion versus frequency constructed from individual data 
points taken at predetermined frequencies. In all tests, the 
excitation was generated in phase by both radial magnetic 
bearings. The excitation orientation was perpendicular to 
the blade mounting plane. An acceleration value greater 
than 1 . 0 g appears at the rotor first-bending critical (355 Hz) 
for each of the four excitation amplitudes. Figure 20 shows 
the corresponding base displacements determined by di- 
viding the acceleration by the frequency squared: 
x = (g X 32.2 ft/s 2 X 12 in./ft X 1000 mil/in.)/co 2 . 


Figure 18. — Spin rig rotor assembly showing location of 
accelerometers on blade and blade base. 


Amplitude, 



Figure 19. — Single-point, blade-base acceleration data 
for bounce excitation (nonrotating shaft). 



Figure 20. — Calculated blade-base displacements. 
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Base acceleration data for excitation current amplitudes 
up to 3 A in figure 19 are roughly constant with respect to 
frequency between 50 and 1000 Hz, and correspondingly 
the resulting base displacement amplitude decreases qua- 
dratically with frequency (fig. 20). The sharp rolloff after 
1 kHz is due to the combined bandwidths of the magnetic 
suspension system components (control computer, filters, 
sensors, actuators). The transfer function (the ratio of 
blade-base acceleration to the magnetic-bearing control 
current) in figures 21 and 22 provides a more thorough 
investigation of how the rotor responded to an excitation 
current amplitude of 1 A. The data were collected by using 
a swept-sine excitation. Figure 21 shows the acceleration 
magnitude response at the blade base. It is clear that rotor 
response is affected by the presence of the blade(s).The 


Blade 

first-bending Blade 



Frequency, Hz 

Figure 21 . — Blade-base acceleration amplitude data 
obtained using swept-sine bounce excitation. 


perception of the authors is that these relatively thick blades 
are heavier and have higher resonances than most blades 
that will be tested. It should also be noted here that these 
tests were performed for one set of controller gains. These 
gains appear to provide the best levitation but are not 
necessarily the ones that will handle all dynamic spin 
conditions. 

The peak at 144 Hz fortuitously matches the blade first- 
bending frequency measured by bench holography (table I) . 
However, in the rig, the expectation is that this mode will 
be shifted and split by two effects. First, the authors expect 
a mistuning effect caused by the mass loading of the 
accelerometer that is mounted near the tip of blade A 
(blade B has no accelerometer). Second, even if the accel- 
erometer were not present, we would expect two system 



Frequency, Hz 


Figure 22. — Corresponding base excitation phase 
angle between swept-sine bounce excitation and 
acceleration measurement. 


TABLE I.— TITANIUM TEST BLADE RESONANCES 


Test 

blade 

Blade mode 

First bending | First torsion | Second bending | First mixed | Third bending | Second mixed | Second torsion 

Resonance, Hz 

A 

144.0 | 667.0 | 879.4 | 2107.5 | 2261.7 | 3753.6 | 4172.9 


Holographic images of resonating blade A 



145.7 



675.9 



900.9 



2140.8 



2301.2 



3827.1 



4258.8 
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Frequency, Hz 


Figure 23. — Blade-tip acceleration magnitude using 
swept-sine excitations of 0.5, 1 .0, and 2.0 A. 



1 00 1 000 


Frequency, Hz 

Figure 24. — Phase of swept-sine magnetic-bearing 
excitation to blade-tip acceleration. 


modes that involve the motions of the two blades coupled 
with the transverse motion of the shaft and the rotation of 
the shaft about its axis. One mode, with the blades moving 
in phase, reacts against shaft translational motion. The 
other, with the blades moving out of phase, reacts against 
shaft rotational motion about its axis. The latter mode can 
only be excited by the magnetic bearings if there is blade 
mistuning. Therefore, there are potentially four modes for 
the mistuned system. However, only two modes are visible 
in figures 2 1 and 22. Because the base-mounted accelerom- 
eter is located on the blade centerline, it is not expected to 
be sensitive to first-torsional motion of the blade; hence, no 
peak is seen near the bench-measured first-torsion fre- 
quency (-670 Hz). A strong peak is seen in figure 21 
around 870 Hz, not far from the bench second-bending 
frequency. The only rotor mode in this plot is that repre- 
sented by the third peak and it does show the largest 
amplitude at 355 Hz. This frequency is the rotor first- 
bending natural frequency as determined by impact testing 
and analysis. The character of the fourth and sixth peaks 



Figure 25. — Blade excitation and free-vibration 
decay at 3000 rpm. 



Figure 26. — Exponential sine fit to decaying free 
vibration shown in fig. 25. 


remains undetermined. Further analysis and measurement 
will be required to identify positively their true nature. 

Figure 22 shows the corresponding phase response from 
excitation to acceleration. The phase lag is substantial, 
confirming that excitation will have to lead the zero-speed 
excitation orientation with respect to the blade row by a 
predetermined amount. This required lead angle will be a 
function of rotor speed and excitation frequency and is 
partly determined by the data in figure 22. 

Figures 23 and 24 present the transfer function curves 
(ratio of blade-tip acceleration to magnetic-bearing control 
current) resulting from a 1 -A, 0-to-peak-amplitude, swept- 
sine excitation perpendicular to the vertical plane that 
contains the blades. Tests at 0.5-, 1.0-, and 2.0-A current 
amplitudes indicate that tip accelerations on the order of 
about 100 g/A are achievable at 0 speed at the blade 
first-bending critical of about 144 Hz and at the first- 
torsional and second-bending-frequency blade modes. The 
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first-torsional mode was excited and captured (albeit not 
using the tilt mode excitation setting) because of the off- 
axis placement of the accelerometer as depicted in figure 1 8 . 

Subsequent to the nonrotating test discussed above, 
blade excitation at 3000 rpm was examined. Figure 25 
shows blade displacements at the first-bending-frequency 
mode during rotation, followed by free-vibration decay. 
The synchronous-rotating-force feature was engaged, thus 
enabling excitation to be maintained at 90° relative to the 
blade row during rotation. These measurements were ob- 
tained with stationary fiber optical displacement transduc- 
ers that sampled each blade once per revolution. Since the 
natural vibration frequency was significantly higher than 
the rotational frequency, the displacement amplitude data 
are aliased about half the rotational frequency. The true 
frequency, however, can be recovered with the use of the 
second transducer, which was installed 30° away from the 
first. The frequency (144 Hz) was obtained from the 
displacement phase angles and the known separation be- 
tween the two transducers. The apparent damping obtained 
from the least-squares fit illustrated in figure 26, therefore, 
must be divided by the ratio of the true natural frequency 
and measured aliased frequency. This ratio was 15 and 
yielded the value of true critical damping of 0.003. One can 
observe from figure 25 that a displacement amplitude of 
70 mils was achieved with this blade sample. 

Concluding Remarks 

The excitation test results of the unbladed and bladed 
rotor clearly demonstrated the usefulness of the Five- Axis, 


Three-Magnetic-Bearing Dynamic Spin Rig as a versatile 
tool for doing detailed analyses of the vibration and damp- 
ing of bladed disks. Factors contributing to this versatility 
and usefulness include (1) the implementation of three 
active magnetic bearings that fully levitate the rotor under 
closed-loop control, thus permitting larger rotor translation 
and tilt displacement and (2) the use of an enhanced control 
code. The code will allow an operator to (a) optimize the 
blade excitation angle that permits maximum resonance 
response from the rotating and nonrotating blades and 
(b) optimize the pivot location during tilt mode excitation 
tests to achieve maximum torsional response. Other code 
refinements include soft touchdown logic, which is used 
for gently delevitating the thrust plate onto the backup 
bearings, and demagnetization logic, which is used for 
demagnetizing the thrust plate. 

A rotor first-bending-frequency mode of 355 Hz was 
discerned from the plot of acceleration versus frequency. 
Blade-tip excitation amplitude (at 0 rpm) of approximately 
100 g/A was achieved at the first-bending critical fre- 
quency mode (-144 Hz) and at the first-torsional and 
second-bending-frequency blade modes. A blade-tip dis- 
placement of 70 mils was achieved at the blade first- 
bending-frequency mode by exciting the blades at a 
forced-excitation phase angle of 90° relative to the vertical 
plane containing the blades while simultaneously rotating 
the shaft at 3000 rpm. 


Glenn Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 7, 2003 
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Appendix A 

Magnetic Suspension Theory 


To calculate the attractive force between the ferromag- 
netic stator and rotor where \i r » 1 , one has to consider first 
the energy W a extant in the air gap between the stator and 
rotor (ref 8). If the field in the air gap is assumed homoge- 
neous as represented in figure 27, the energy W a stored in 
the volume V a , has the form 


W a =-B A H a V a =-B A H a A a 2s 

a 2 a a 2 a ^cl a 


(Al) 


It is evident from equation (A3) that the magnetic force 
has a quadratic dependence on current and an inverse 
quadratic dependence on the air gap. The net controlling 
force (due to an opposing pair of identical electromagnets) 
acting on a rotor can thus be represented as (ref 5) 


F=z 





(A4) 


A change in the field energy in the air gap with respect to 
rotor position produces an electromagnetic force. This 
force is derived in equation (A2) wherein the relation B = 

p^p^i/is used and the magnetic flux F^ a A a is considered 
invariant for infinitesimal displacement ds, hence, 


, „ , B 2 A A a 

f„= — -=b a h a a = — ^ — 

ds a a (i 0 


(A2) 


where q and z 2 are the current in the opposing coils; sq and 
s 2 are the gap distances between the rotor and each oppos- 
ing pole face. 

The squared terms in equation (A4) are undesirable from 
a control standpoint and are eliminated by using a lineariz- 
ing technique that incorporates a bias current and a control 
current. By replacing i x and z 2 in equation (A4) with 
(i h + Q and (i b - z' c ), respectively, and sq and s 2 with 
(^q - s de ) and (sq + s de ), respectively, the force exerted by 
the two actuators takes the form 


Inserting equation (B6) in equation (A2) produces the force 
equation: 


fa 


f ni 


,2 


,2 


\2sj 

where z = (p i Q n 2 A a )/4. 


1 2 A l l /A 

4 / 



Figure 27. — Rotor-stator segment. 


F=z 


{jb +z c ) (jb i c ) 

2 2 
( s 0~ s de) (^0 +s de ) 


(A5) 


where i b and i c are the bias current and the control current, 
respectively; s 0 is the nominal gap distance, and s de is the 
deviation from the nominal value. After making the appro- 
priate algebraic manipulation and taking the requisite par- 
tial derivative, the linearized force current, and position 
are shown to have the approximate relationship 

Fi^Sde+Xfic (A6) 

for i c « i b and s de « s 0 and where X s is the position 
stiffness and X t is the current stiffness. In the case of 
proportional-derivative feedback control when an excita- 
tion signal is used, z c is replaced with ~{^K p s de +K d s de }+i ex 

where K p and K d are the proportional control gain and 
derivative control gain, respectively, and i ex is the current 
used to produce a forced excitation of the shaft. The 
equation of motion for the shaft displacement s de is 
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^ex m eq $de *^i s de ^s) s de (^ 7 ) 

where m is the rigid rotor equivalent lumped mass and F ex 
= Further algebraic simplification produces an ex- 
pression of the form 


^ex m eq $de c eq s de ^eq s de 


+K 


(A8) 


Therefore, 


m 0 


eq s de ^ eq s de c eq s de+^ex 




(A9) 


where c eq = X )K d and k = TkK - 7k, . 

Equation (A9) is at the heart of magnetic-bearing control 
and is the crux of the five-axis-rig control code (ref 6 and 
its latest version). 
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Appendix B 
Magnetic Circuit 


In the case of a magnetic bearing where iron is used in the 
stator and rotor (ref. 8), electromagnets cause a flux <E> to 
circulate in a magnetic loop (see fig. 28 for an approxima- 
tion to the flux path through the stator and rotor). One 
assumes that the flux is entirely confined to the iron (no flux 
leakage) except in the air gap. Further assumed is that the 
flux runs entirely within the magnetic loop in which the 
iron cross section Aj e is constant along the entire loop and 
is equal to the cross section A a in the air gap. Thus, the 
relation can be stated as 

®=B Afl Afe= B Aa A a (Bl) 

where A^ e = A a and 


(B— ) 

For a homogeneous field within the magnetic loop, 
consider the mean length L of the magnetic displacement 
and an air gap of length 2s. If the iron has constant 
permeability, the equation for the magnetic circuit depicted 
in figure 28 can be written as 

j>Hds=l fe H Afe +2sH Aa =ni (B3) 

The term ni in equation (B3) is usually referred to as the 
“magnetomotive force” (mmf). It is evident from equation 
(B2) that the flux density B in the iron and the air gap is 

identical. The field intensities H Af e and H Aa inequation 

(B3) can be replaced by using the relation B = 
hence, 



B 


B 


l fe 


+2 s — —ni 
M'OM'r M'O 


(B4) 


Explicitly solving for B yields the expression 


_ m 

B =H~, , 

—+2s 

V^r ) 


(B5) 


For iron in which p, » 1, equation (B5) simplifies to 


n m 

S=n°- 


(B6) 
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